A "quasi-complete" mechanical model for a double torsion pendulum 
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We present a dynamical model for the double torsion pendulum nicknamed 'PETER', where one 
torsion pendulum hangs in cascade, but off-axis, from the other. The dynamics of interest in these 
devices lies around the torsional resonance, that is at very low frequencies (mHz). However, we find 

^•f^ that, in order to properly describe the forced motion of the pendulums, also other modes must be 

considered, namely swinging and bouncing oscillations of the two suspended masses, that resonate 

f"-) at higher frequencies (Hz). 

£Xj Although the system has obviously 6+6 Degrees of Freedom, we find that 8 are sufficient for an 

, | accurate description of the observed motion. This model produces reliable estimates of the response 

to generic external disturbances and actuating forces or torques. In particular, we compute the 
effect of seismic floor motion ('tilt' noise) on the low frequency part of the signal spectra and show 
that it properly accounts for most of the measured low frequency noise. 

PACS numbers: 04.80.-y, 04.80.Nn,95.55.Ym, 

cr 

jll I. INTRODUCTION 

Free Fall, i.e. motion in absence of external forces, is a key ingredient of many present [1] and future |3] space 
i-H missions. R is particularly important for proposed LISA-like missions, that aim to detect gravitational waves by 
measuring the change of distance between two Test Masses in geodesic motion. Indeed the precursor LISA-Pathfindcr 
m e] win soon be launched to demonstrate many of the key technologies that LISA is based on, including the 
effective realization of a free falling Test Mass (TM), reducing the residual disturbances to the level of lOfN/^/Hz at 
lmHz. Such demanding requirements need ground based facilities for preliminary free fall tests. This has prompted 
an extraordinary effort into reducing the effects of gravity and other local forces. The ideal tool for these tests is, 
traditionally, the torsion pendulum, where gravity is balanced by a supporting fibre and rotational motion is virtually 
uncostrained around it. Indeed, the group at the University of Trento has extensively tested residual forces on two 
different apparata: a torsion pendulum where the TM moves freely in rotation [B] and one where the "soft" motion is 
approximatively a linear translation T . The opportunity to experiment with a TM being in quasi free fall over more 
than one degree of freedom (DoF) is quite desirable, as it better represents geodesic behavior and it allows investigation 
on possible spurious couplings between motion on different DoFs. This is crucial for the drag-free operation of the Test 
Masses [5] on LISA-Pathfinder (and of any of the possible realizations of the LISA concept): in such space mission, 
the TMs will move freely in the sensitive direction, but will be under feedback control on the remaining 5 directions. 
An instrument with two soft DoF, like a double torsion pendulum can observe the free motion of one DoF when the 
other is actively controlled and can provide useful information on the amount of cross-talk that the control electronics 
might feed on the "free" mode. 

To this purpose, we have developed an instrument where one torsion pendulum hangs, off-axis, from another. This 
results in the TM being almost free (or having soft restoring forces) both on a rotational and on a translational degree 
of freedom. The dynamics of this system is quite more complex than that of a single DoF torsion pendulum and its 
description requires a careful modelling of all its (2 x 6, in principle) DoF, because some of the external disturbances 
(tilt motion, for example) couple to torsion via higher frequency modes, that need be accounted for in the model. 
We present an analytical mechanical model of the double torsion pendulum, taking as reference the instrument built 
and operated in the INFN lab in Firenze (I), nicknamed PETER (PEndolo Translazionale E Rotazionale, namely 
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Translational and Rotational PEndulum) [91 110) . 

In this instrument, a cubic TM (see the left panel in figure [T]) is suspended through a double stage system with two 
torsion fibres in cascade: the lower fibre allows this mass to be almost free in rotation around its symmetry vertical 
axis. The fibre hangs from the tip of one arm of a crossbar that is, in turn, suspended with an upper fibre to the 
fixed frame, represented by the vacuum enclosure. The torsion of this upper fibre allows almost free motion of the 
suspended system, including the test mass, along an arc of a 30 cm diameter circumference: for small torsion angles, 
this can be considered for all practical purposes a translational motion. Dummy loads hang from the three other arms 
of the crossbar. The TM is enclosed in the Gravitational Reference Sensor (GRS) a hollow metal box padded 
with electrodes that permit to monitor the motion of the TM along all its translational and rotational DoF. 

The TM is a hollow Al cube with a 46 mm edge plus a shaft (81.5 mm long) that connects to the fibre and 
provides electical insulation. The cubic TM (but not its shaft), the GRS and the readout electronics, well match 
the set-up of the LISA-Pathfinder flight model geometry. The apparatus is also equipped with additonal readouts: 
an autocollimator and an Optical Read Out (ORO) system [T2] that provide independent measurements of the test 
mass along the 2 soft DoF. The sensitivity goal for this apparatus, when limited by the mechanical thermal noise and 
the readout noise, is better than 10~ 13 m s _2 /Hz 1 ^ 2 around lmHz (on each DoF), namely, only 1 order of magnitude 
worse than the LISA-Pathfinder goal along the sensitivity axis. A more detailed description of the apparatus can be 
found in [IH]. 

Aim of this paper is to develop an analytical model of the whole system by which its main features can be evaluated 
and compared with preliminary data runs. In the title the quotation marks around "quasi-complete" stand to point 
out that the model goes as far as possible in the comparison with the actual experiment. Although the operation 
range is in the mHz band around the torsional resonances, the model includes the treatment of the swinging pendulum 
and bouncing resonances, that take place at much higher frequencies (~Hz). 

Such a detailed description may appear superfluous: however, we have carried out the analysis of the mid-high 
frequency range for an overall validation of the model and with the purpose of analysing possible sources of external 
disturbances. In particular, in view of the complex structure of PETER, the tilt noise at low frequencies can be 
properly described only if the coupled swinging pendulum motions are taken into account. 

The plan of the paper is as follows: in sect. |ll]we lay out the mechanical model, with special attention to defining 
reference frames that are suitable to describe motion of the two payloads (the crossbar and the Test Mass) in the 
appropriate limit of small oscillations. The dynamical properties of the double pendulum, i.e. angular velocities and 



moments of inertia are then introduced and, from those, the Lagrangian of the system is derived. Sect. Ill analyizes 
the free motion of the double pendulum and derives the normal modes and their resonant frequencies, while in sect. 
IV[ introducing the generalized forces on the system, we discuss how the influence of external disturbances on the 
system can affect the pendulum output. 



II. EIGHT DEGREES OF FREEDOM MECHANICAL MODEL 



In principle, in order to describe a double pendulum of the PETER kind, with two "payloads" suspended in cascade 
through two torsion fibres, we should consider 6 + 6 DoF, to allow each payload to translate and rotate in every 
possible way. This, in the constrained motion of a pendulum, maps into describing the torsion, 2 pendulum motions 
(in 2 orthogonal directions), bouncing of the fibre length and 2 rocking (rotations around a horizontal axis passing 
through the suspension point) motions. However, in all data gathered so far, the rocking modes, expected at the 
frequencies of 0.55Hz and 3.4Hz, were not detected: if they exist, their amplitude is well below the noise level. 
Therefore, in our analysis, we neglected the rocking motion of both the crossbar and the TM, assuming that they 
rigidly move with the respective fibres that, therefore do not bend at the suspension point. While their inclusion 
would lead to more cumbersome expressions for the equations of motion, as a matter of principle they can be included 
without problems. 

The 'quasi-complete' model for PETER describes then the torsional and pendulum oscillations of the crossbar and 
the TM, including also their 'bouncing' motions. Therefore we deal with a 4+4 = 8 DoF model. 



A. Reference frames and configuration variables 

We need to describe the position and rotation of two rigid bodies (the crossbar and the TM) with respect to an 
inertial frame. However, we find it useful to introduce also two moving (non-inertial) frames and exploit the Lagrangian 
formalism to combine generalized coordinates for a simpler description of the dynamics. The two additional frames 
are comoving with the two load masses and will be labeled with the "a" (above) and "b" (below) subscripts. In 
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FIG. 1: Left panel. Schematic representation of PETER with the three triads of unit vectors. Right panel. The 
coordinate systems described in the text: u Xl u y , u z represent the inertial frame; P is the position of the centre of 
mass of the crossbar, Ui, Ui and U3 is one of the moving coordinate systems; 9, 77 and ip the three angles to identify 

the orientation of the frames. 



analogy with the treatment of the 1 DoF case |13j . we have made a choice of the angles suitable to describe small 
oscillations around the stable equilibrium of the whole system [201 HI] • 

We set the origin of the inertial frame in the suspension point of the upper fibre. The orientation of this frame is as 
follows (see the left panel of fig[l]): u z along the local gravity acceleration — g, u x and u y along the crossbar arms at 
the equilibrium position (when no torsion is present in the upper fibre). The first non inertial frame, comoving with 
the crossbar (the "a" frame, Ui a , u 2 . a ^ u 3 a ) has its origin at the end of the crossbar arm from which the lower fibre 
is suspended; Ui^ a is parallel to this arm and directed towards the centre, U\^ a is parallel to the other arm. Since we 
assume no rocking motion, u^ ta is at all times parallel to the upper fibre direction. 

The second moving frame (the "6" frame, whose origin coincides with the center of mass of the TM) describes 
the orientation of the TM with respect to the "a" frame (i.e. to the crossbar). The unit vectors «i 1*2. 6, 1*3,6 are 
directed along the principal axes of inertia of the TM, and therefore u^^ is parallel to the lower fibre direction. In 
order to bring the "a" frame with its axes parallel to those of the inertial frame, we introduce the following sets of 
rotations (see the right panel of fig[l]): first, by an angle 9 a around the x-axis; then by an angle r\ a around the new 
y'-axis; finally, by an angle ip a around the new z'-axis, now the same as u 3 . The resulting orthogonal rotation matrix 
R a is 
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cos rj a cos (fa sin r\ a sin 9 a cos <p a + cos 9 a sin <^> a sin 9 a sin </3 a — sin r\ a cos a cos ip a \ 
R a = I - cos ?] a sin <^ a cos # a cos v? a - sin r\ a sin # a sin v? a sin r\ a cos # a sin <p a + sin 6> a cos (p a \ (1) 
sinrya -cos?7 a sin# a cosr? a cos# a y 

Analogously, a second rotation by the angles 9 bl ?7&, <fb, denning a matrix _Rfc, is used for the transformation a — > b 
between the two moving frames. 

We now introduce an additional convention: the superscripts ' and " indicate that the components of a certain 
vector are expressed in the a or & frame, respectively. An absence of superscripts indicates that the vector is defined 
in the inertial frame. In the inertial frame, the unit vectors t*i i0 , u 2ja , u 3 a are: 

ui >a =^(1,0,0) T ; U2,„=J£(0 > 1,0) T ; u 3 , a = fl£(0, 0, 1) T . 

The unit vectors at rest with the TM are given, in the "a" frame, by 

u' li6 = JI? , (l,0,0) T ) u' 2tb = RU0,l,0f, u' 3tb = Rl(0,0,l) T 

while the same vectors are expressed in the inertial frame as: 

Ui i, = RaUi b with £ = 1,2, 3. 

We now need to express the position of the center of mass, r g «, of the crossbar system and r g _ b of the TM in the 



inertial frame (see fig. [2]). 
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•G b 



FIG. 2: Schematic representation of the heights of the centers of mass of the double pendulum. G a and G b are the 
centers of mass of crossbar+counterweights and of the TM, respectively. 

We define Z as the component along u 3 ^ a of the distance between the fibres suspension points, Z a as the distance, 
again along tt3 >a , between the the lower fibre suspension point and the center of mass of the crossbar system. S a (t) 
describes the elongation of the upper fibre. We get 



9,a 



-(Z + Z a + S a (t))u 3 



(2) 
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where y a is the horizontal displacement of the center of mass (with respect to the line of the upper fibre) . It can be 
easily computed by assuming no misalignments: in this case the center of mass of the whole system (at the equilibrium 
position) has null y-coordinate, therefore 

,m b 

Va = d — 
m a 

where d is the length of a crossbar arm, m a is the total mass of the crossbar system and mt is the mass of the TM. 

Finally, if Z b is the distance between the lower fibre suspension point and the centre of mass of the TM and 5 b (t) 
is the elongation of the lower fibre, for the Test Mass we obtain 

Tg.b = -{Z + S a (t))u 3ia - du 2 ,a ~ (Z b + S b (t))u 3ib (3) 

B. Angular velocities 

Using the Poisson formula for the time derivative of a unit vector, one can readily derive the angular velocities of 
the crossbar system and of the TM in the fixed frame: 

3 3 

1 ^ du lM 1 ^ dui, b 

"a = 2 XJ «,,„ X — , 0* = - £ U t , X — 

i— 1 i— 1 

In the moving frames, they are respectively transformed to: 

Wi = Ra "a, ^fc = RbRa ' V b . 

C. Inertia tensors 

We denote with I a the inertia matrix for the crossbar element, containing three counterweights and their rigid 
connection to the crossbar, and with I b that for the TM element, which includes the TM as well as the shaft connecting 
it to fiber b. Both inertia matrices are calculated with respect to their centers of mass and in the references frames 
(iti.a, U2, a , u 3,a) and («i,6 , 1*2,6, u 3,b) , respectively The plane x = is a symmetry plane for the crossbar system, 
therefore 7f 2 an d ^13 are expected to be zero or, at least, quite smaller than the others coefficients. Instead, the inertia 
matrix of the TM is expected to be diagonal and with two terms equal to each other. 

flii \ //fi 

la = I 2 a 2 /fa , /t= 4 

V J 23 J V 7 3 b 3 

D. Lagrangian and generalised coordinates for the system 

A rigid body with inertia tensor Ijk, j, k = 1, 2, 3 and angular velocity ujk has rotational kinetic energy 



K rot = \ ^2 I kUjU k = \u ■ I ■ 



2 

jk 

The translational kinetic energy is instead obtained from the velocity of the centre of mass: by using the above defined 
coordinates, we therefore obtain: 

K — — I ' I 2 K — — I ' I 2 

^transl.a — '^' m a\^g,a\ 5 ^-transl.b — l^^^^Q^ 

and 

K ro tM = ' la ■ W a , Krot.b = ~j^ b ' h ' U b ■ 
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The total potential energy is 



U = m a gz gta + m b gz g .b + \^k a Lp 2 a + ^k b Lpl + -K e<a S 2 a + ^K e<b S b , 

where z g ^ a and z g ^ b are respectively the heights of the centers of mass of crossbar and TM, k a and k b are the torsional 
constants of the two fibres and K e ^ a and K e . b are the bouncing spring constants, responsible for the fibre elongation 
[T3]. The fibres are considered here as ideal elastic wires, as discussed in sect IV B| Therefore, we neglect bending 
contributions to the potential energy. 

The total Lagrangian of the system is then 

M?! Q) — KtransLa + K trans i b + K rot a + K rot b — U, 

where with q we denote the set of configurational coordinates 6 a , T] a , (p a , S a , b , r) b , ip b) S b . This Lagrangian is only 
apparently decoupled in its variables: indeed, the coupling is hidden in the definition of r 9ib (eqj3|. We now exploit 
the limit of motions of small amplitude, and construct a quadratic Lagrangian from which we calculate the equations 
of motion (reported in Appendix). 

In the limit of small oscillations, the observable quantities (position and orientation of the TM as measured by the 
GRS) are related to the Lagrangian coordinates by 

' -r} a (Z + Z b + h c ) - j] b (h c + Z b ) + ip a d\ 

6 a (Z + Z b + h c ) + 6 b {h c + Z b ) | ; (4) 
S a - $b — d a d 



e\ fe a + e b \ 

rj ~ ri + 776 

<Pj \tPa + <PbJ 

where h c is the distance between the cube geometrical center and the center of mass of the TM torsion member (i.e. 
G b in Figure [2| , which are offset due to the mass of the shaft connecting the TM to the torsion fibre. 
We remark here that <p a is an additional observable quantity, as it can be monitored separately, e.g. via an au- 
tocollimator measuring rotation of the crossbar. Moreover, it is a good observable, because it relates in a simple, 
straightforward way to any linear force acting on the TM. 



III. NORMAL MODES AND FREE MOTION 



To validate the model we now compare its predictions with real data from the PETER apparatus [22]. Lengths, 
masses and inertia moments are measured or computed with good accuracy. Springs constants are derived from 
measurements of resonant frequencies of each oscillator. However, when analyzing a 2 DoF system, we observe the 
resonances of the coupled modes, rather than those of the physical oscillators. For this reason, we used preliminary 
data when the crossbar was clamped and we effectively only had the isolated "b" pendulum, with the TM moving 
along four DoF (2 swinging modes, 1 torsional mode and 1 bouncing mode). We measured the resonant frequencies: 

^torsional — 2.2 HlHz, V b ouncing — 8.82 Hz . 

From these measured values, we derived (see [13] for details) the torsional and bouncing spring constants of the "6" 
fibre, reported in table 1. We then estimated the spring constants of the "a" pendulum, by applying the following 
scaling relations, involving the radii (r a ,r b ) and lengths (/ a ,/t) of the fibres,: 




The lengths of the fibres are re-defined by taking into account the static longitudinal deformation (a few mm in both 
cases) due to their respective loads. The torsion constant of the 'a" pendulum was measured on an independent, 
dedicated test apparatus, yielding a value in excellent agreement with the above determination via scaling. All 
physical and geometrical parameters of the double pendulum are summarized in Table [I] measurement errors amount 
to a few percent. By using these values we compute the frequencies of the normal modes, shown in the second column 
of Table |IT] A least square fit, constrained within the experimental error bars, then adjusts the mechanical parameters 
in order to best match the observed normal mode frequencies, reported in column 3 of the same Table. 
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Masses [kg] 




m a =1.2 


m 6 =0.11 


Torsional constants [kg 


2 —21 
m s J 


k a = 1.8- 10~ 6 


fc 6 = 7.1 • 10" 9 


Bouncing constants [kg 




H e _ = 4804 


K e ,b = 300 


Lengths [cm] 




d=15; h c =3.4; 2=87; 


2 b =76; 2 a =43 



TABLE I: Numerical vaues of the mechanical parameters of 
the PETER double pendulum. Masses and lengths are measured, inertia moments and spring constants are derived from 
measured quantities. Experimental errors are of the order of 2-3 % for most values, and 10% for the torsion and bouncing 
constants. These quantities are used as input parameters for the model, yielding the eigenfrequencies listed in table [III 



Normal modes frequencies 


mode # 


calculated 


measured 


unit 


V\ 


1.3 ±0.1 


1.331 


mHz 


V2 


2.2 ±0.2 


2.117 


mHz 


V3 


0.41 ±0.01 


0.406 


Hz 


V± 


0.42 ±0.01 


0.4065 


Hz 


V5 


0.59 ±0.01 


0.58815 


Hz 


V6 


0.62 ±0.01 


0.6170 


Hz 


v-r 


8.05 ±0.6 


7.926 


Hz 




10.9 ±0.8 


10.393 


Hz 



TABLE II: Calculated and measured normal modes frequencies. The error on the calculated values are obtained by 
a Monte Carlo variation of the input parameters within their experimental error, given in table IT) 



The two pairs ^34 (virtually degenerate) and v^^ are respectively associated to the swinging motion of the cross- 
bar±counterweights system and of the TM. Indeed, if each pendulum behaved as a decoupled physical pendulum, wc 
would get 

/ m a g(Z + Z a ) . ..„„ 

V It + m b Zf; 



There are two main reasons for the line splitting v 5 ,v§, associated with the swinging of the TM: the first is that 
the centers of mass of both the crossbar±counterweight system and of the TM are not aligned with the upper fibre 
(cfr. the case of 1 DoF, [12]); the second is that the inertia matrix of the crossbar±counterweight system is not 
diagonal. In other words, the "b" pendulum couples differently to the crossbar system, depending on whether the 
lower pendulum swings parallel or normal to it. We then solved the equations of free motion (namely the homogeneous 
system of equations associated to (A2|) with many arbitrary but 'reasonable' initial conditions. From the solutions 
we have obtained the observable coordinates by using Eqs. (IID). The solutions were then sampled at 0.1 s, just as 
in the actual experiment, and from these the spectra of the simulated data stream were computed. A white noise of 
dnm/y^Hz) was added to the model spectra, to simulate the readout noise. Note that damping is not considered 
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FIG. 3: Amplitude spectral densities of observable displacements X, Y, Z. Left: in the frequency region around the 
torsional resonances. Rigth: in the intermediate frequency range, where the swinging resonances are observed. 
Upper plots show the model predictions, while the lower plots display the experimental data. Solid line (black 

online): X, dashed line (blue online: Y, dotted lines (red online): Z. Only the motion along X is activated at the 

torsion eigcnfrcqucncies. 



in the model; therefore, the width of the resonances is only determined by the length of the time series, that roughly 
corresponds to 12 hours of data. Finally, we compared these spectra with those of preliminary experimental runs, as 
shown in fig. ([3]). The experimental spectra have a much larger wide band noise that is partly due to tilt effects and 
partly under investigation. Similar comparisons have been carried out for the rotational observables (tp, r\ 9) that are 
measured by the GRS. 

We note that the model well replicates most of the features of the actual experiment, including the non trivial 
prediction about which resonances appear in a given observable channel. We verify, for example, that the swinging 
resonances of the lower fibre, at 0.6 Hz in fig[3| are split in a doublet (v 5 and i/g), as stated above. The X(and 77, 
not shown) channels only see the vq mode, and that is correctly predicted by the simulation. However, discrepancies 
remain: there are modes that show up in some channels in the data, while not predicted as, e. g., v§ in the Y and Z 
spectra. We recall however that the comparison is not completely fair: simulated data only predict the free motion 
of the double pendulum, while the instrument is certainly driven by external disturbances. Besides, unavoidable 
asymmetries in the assembly of the system (e.g., an imperfect match of weigth on the four arms of the crossbar, that 
can lead to the crossbar laying at rest in a non-horizontal plane) are not accounted for in the model, and can easily 
lead to the appearance of these modes 

In the following section we discuss how to include driving effects, through the formalism of Lagrangian generalized 
forces. We apply this method, as an example, to one of the most obvious and easily measurable of these external noise 
sources, i. e. floor tilt, and show how this forcing term, even if only relevant at very low frequencies, can produce 
excitation of higher modes. 
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In order to further validate the model, we have undertaken another comparison: the TM was moved away from 
equilibrium by a large kick along the X axis by the electrostatic actuation, and then released. We obtain a good 
approximation of free motion, as possible disturbances are only relevant at smaller amplitudes. The same initial 
condition were then used for the model, and the resulting motion calculated. Figure [4] shows the comparison of real 
and simulated data for the tp channel of the GRS: while the match of the low frequency behavior (a sum of the two 
normal mode oscillations at v\ and z^) is excellent, the measured data show a modulation at the swinging frequencies 
that the model does not replicate. However, figure [5] shows a zoom on a smaller stretch of time of the X observable, 
showing a good agreement also for the motion at the swinging frequencies. 




500 1000 1500 2000 2500 3000 
time [s] 

FIG. 4: Measured (dashed-blue) and simulated (continuous-black) cp evolution of a large amplitude, virtually free 
motion: low frequency behaviour of the model well replicates the experimental data. The slight mismatch in the 
amplitudes of vibration can easily be due to imperfect matching of the gain calibration coefficients (±2%) or 

experimental errors on the model parameters. 




time [s] 

FIG. 5: Measured (dashed-blue) and computed (continuous-black) X evolution of of a large amplitude, virtually free 
motion. A fraction of a torsion period is shown, in order to better display the motion at the higher, swinging 

frequencies. 



IV. GENERALIZED FORCES 



PETER will be not, in general, an isolated system. External disturbances such as the seismic vibration of the 
laboratory, parasitic electrical or magnetic fields, mechanical dissipations etc. can be modelled as time-dependent 
forces and torques applied to the crossbar and/or the TM. Moreover, by changing the voltage on the electrodes one 
can induce on the TM a deterministic force and/or momentum (actuation). Actuations can be periodic functions, 
but it is also possible to induce a 'viscous damping' in order to reduce the amplitude of the oscillations of the TM. 



A. Small external disturbances 



In the most general case, two forces (-F , Fb) and two torques (M a , Mb) with arbitrary intensity and direction must 
be considered. Momenta are calculated with respect to the centre of mass of the crossbar+counterweights system and 
the TM. In the limit of small displacements, we can perform a series expansion and keep the lowest order non-trivial 
terms. 

In the general case of a non-isolated system of M particles with N DoF (N < 3M if there are constraints, in our 
case, simply, N = 2), we can write the Euler-Lagrange equations in the presence of generalized forces: 

d dC dC _ „ ^ „ dr, , 

*«^-«- = ^i Qh = y\* , rjr A - h=i,...,N 

dtdq h dq h f-J dq h 

where Fj is the force applied on the j-th particle. If the system consists of a set of M' rigid bodies, the generalized 
forces are written as: 



i=l 



(5) 



where the Ft < and Mi are respectively the resulting external force and momentum and the tyi are implicitly given 
by 

We apply eq.|5| to our case [M 1 = 2), with two generic forces F a , Fb and torques M a , Mb applied on the TM and 
the crossbar. We report in Appendix the explicit expressions for the Qh, linearized in the configuration variables. 



B. Seismic effects 



Generalized forces can be used to perform an analysis of seismic effects on the motion of PETER. With respect to 
an ideal inertial frame, the floor of the laboratory is subject to small translations and rotations ("tilt" noise). With a 
simple model we calculate how the observables of PETER are affected by the motion of the ground. We assume, for 
simplicity, the vacuum chamber where PETER is housed to be a rigid body (this is certainly true at low frequencies) 
and we neglect the flexural stiffness of the fibres [2"5] . 

The tilt of the ground with respect to a "fixed" frame, is describer by an angular velocity that we express, for small 
angles, as u) t = (0t(t) , f)t(t) , <pt(t)) . We focus here on rotations, neglecting for the moment the translation of the 
ground, as well as the negligible rotation on the horizontal plane, <pt(t). Indeed, (9 t (t) , rj t (t)) are directly measured 
by a tilt-meter and can readily provide a test for the model. The whole (rigid) structure rotates around its base with 
angular velocity u)t and, consequently, the upper fibre suspension point (that is at a distance h from the ground) is 
affected by a linear acceleration r t = h{f] tl —6 t ,0). 

We have derived the Lagrangian of the double pendulum in an inertial frame. Since gravitational potential energy 
depends on the elevations of Test Mass and crossbar, the coordinates 9 a and rj a must be corrected, by adding the 
corresponding tilt angles, as 6 a + 9 t and r\ a + rj t , while all other coordinates remain unchanged. After this change in 
the Lagrangian, we must add to the resulting equations of motion the external (apparent) forces due to acceleration 
of the origin and to angular velocity. In the resulting, well known expression for the force on a generic mass element 
dm, of coordinate r, in an accelerated system: 

dF app = —dm(r't + w t x (u) t x r) + 2u t x r + w t x r) (6) 
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we can safely neglect, for small u>t and small oscillations, both the centrifugal and the Coriolis term that represent 
second order corrections. The only relevant terms are then the first and the last one. 

We now compute the total moment of the apparent forces with respect to the center of mass r g of a rigid body, by 
integrating over its volume. With the substitution r* = r — r g , and making repeated use of the property J r* dm = 0, 
we obtain: 

M app = J r* x dF apPtl = - J r* x (r t + 0J t X r)dm. 

As the first term is zero, 

M app = - J r* x (o> t x r g )dm - J r* x (u> t x r*)dm = -I ■ w t 

where I is the inertia matrix of the rigid body. 

In this way, we have derived the apparent forces and torques acting on the Test Mass and crossbar: 

F apP: a = -m a (r t +dj t x r 9:Q ); F apPib = -m b (r t + d) t X r S)6 ); 
M apP:a = —I a ■ u>t; M app . b = —I b ■ d> f . 

Inserting these expressions into eq.|5]), we can now explicitly compute the effect of tilt motion on the torsion pendulums 
observables. 

We report in Figure [6] the transfer functions from tilt angles (6t, rjt,<pt) to the X , ip and tp a coordinates that monitor 
the "soft" DoF and are therefore of special interest. We observe that, when the frequency approaches zero (below the 
torsional resonances), X — > —r\ t (Z + Z b + h c ) and, analogously (not shown), Y — > Qt (Z + Z b + h c ). All other transfer 
functions tend to zero at low frequency, and therefore do not affect the observables in the range of interest. 

Using the above described transfer functions, too cumbersome to be written here, it is possible to calculate the 
spectra of the observables (in particular, X), as expected from tilt measurements. This allows us to assess the role of 
tilt in the low-frequency (below torsional resonances) noise, and its predicted effect on the X variable. Tilt is indeed 
responsible for the observed high level of off-resonance noise in the measured spectrum 

Incidentally, we note that the torsional observable ip responds to tilt motion r\ t at the swinging modes (y^ -j- Uq): 
this is, probably, the reason why these resonances are observed in the (p data, while not predicted by the solution of 
free motion. 

An analogous analysis can be carried on to evaluate the additional effect of translational floor motion. This is not 
as interesting because these transfer functions would all vanish toward zero frequency, when the suspension point 
rigidly translates together with the entire enclosure and, in particular, with the GRS. Besides, as there is no effective 
way to measure such motion, calculation of these transfer functions is of little use. 

C. Low frequency behaviour: the 2 soft DoF 

In many practical cases, we are mainly interested in the low frequencies dynamics of the system, i.e. in the region 
around its two soft, torsional resonances, namely between 1 and 3 mHz. It can then be useful to simplify our model 
in the limit where all swinging and bouncing motions can be neglected: we are left with the two torsional motions 
described by the two angles tp a and (fb, and governed by the equations: 

'Pa {I33 + d Va [m a + m b ) + /f 3 ) + I^ipb + k a ip a = Q Via (7) 

hsVa + ^33^6 + h<Pb = Q v ,b 

with Q Va = M a<z + dF b . x + M b<z , and Q Vb = M b>z . 

It might be useful to re-express the equations in terms of the observables measured by the GRS (and defined in the In- 
ertial System). In the limit considered, where the crossbar and TM move rigidly together (except for torsions), we have 

Xgrs = fad p = (p a + tp b (8) 
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FIG. 6: Transfer functions for the observables X (left panel), (p (right panel) and (p a (bottom panel) relative to the 

ground tilt angles 6 t (solid lines), rj t (dashed), ip t (dotted). 



The corresponding quadratic Lagrangian is 

C T = -I c <j> 2 a + ^I t if 2 - ^Kvl - \h{Va - f) 2 , (9) 

where I t = /| 3 and I c = I| 3 + dy a (m a + nib)- 

The (forced) Euler-Lagrange equations obtained from Ct are 

Ic<Pa + K<Pa ~ h{V ~ <Pa) = M a, (10) 
I t <p + k b (<p-<p a )=M b . (11) 

This simplified model is adequate to describe most of behaviour of the double pendulum at frequencies near or below 
the torsion resonances. However, its predictions might result inaccurate in instances where disturbances feed into the 
low frequency part of the spectrum through coupling with higher modes, like in the example of the previous section. 



V. CONCLUSIONS 



We have described an 8 DoF Lagrangian model that provides a suitable account for the motion of the double torsion 
pendulum "PETER" . The model fully describes its free dynamics and its response to external disturbances, and can 
accurately predict the torsional, swinging pendulum and bouncing resonances. The model is otherwise 'necessary' to 
obtain a correct account of external forces and torques acting on the system. In particular, it gives a good account of 
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the tilt noise, which feeds into the low frequency motion of the double pendulum via a non-trivial coupling among the 
various DoF. The number and location of resonance peaks are correctly predicted and are used for a first validation 
of the model with preliminary data. The transfer functions concerning seismic noise components allow us to evaluate 
the response of the observables, with particular emphasis on the 'soft' translational motion. Dissipation effects were 
not considered here, in order to keep the equations manageable. However, it is possible to extend the Lagrangian 
formalism, via the Rayleigh dissipation function, to account for frictional forces. This is actually straightforward in 



the low-frequency, two DoF limit (cfr. sect. IV C I, where dissipative effects are taken into account in order to predict, 



via the fluctuation-dissipation theorem, the fundamental limits of sensitivity of the system |22j . 

Appendix A: Equations of motion in Lagrangian coordinates. 

We define the following quantities, that are, in essence, modified moments of inertia (the first six are, e.g., referred 
to the upper suspension point). The last term is proportional to the z coordinate of its center of mass. 





= Ifl ' 


fm„ [(z + z a ) 2 + y 2 a )}, 


%b,l 


= I b n - 


- m b [(Z + Z b ) 2 + d 2 )} , 


Ax, 2 


— 1 22 


^m a (Z + Z a ) 2 , 


%b,2 




-m b {Z + Z b ) 2 , 


Za,3 


= ^33 " 




%b,3 


— 1 33 


- m b d 2 , 


%b,4 


= Ill - 


"rn b Z^, 


I a ,4 


= I 23 ~ 


t- m b d(Z a - Zb) , 


%b,5 


— 1 11 


- m b Z b (Z + Z b ) , 



ci = m b (Z + Z b ) + m a (Z + Z a ) . 

By assuming i| 2 = I 12 = -^13 = 0; the equations of motion for the double pendulum are: 

(la.j +l bA )8 a + l b , 5 9 b + m b dS b + 

9 c l@a + ' m b9Z b Q b — Qg^ai 
(Ia,2 + Zb,2)ija + ZaA&a + ^b,5Vb + 9Cl1]a + m b gZ b 7] b = Q V m, 

(I a , 3 + T b ^)ip a + I% 3 tp b +l aA rj a - m b dZ b rj b + k a ip a = Q v>a , 

(m a + m b ) S a + m b S b + <S a K e , a = Qs,a, 

l b , 5 9 a + T bA 9 b + m b gZ b (9 a + 6 b ) = Q e , b , 

l-b.^a +Zb,4Vb - m b dZ b Cp a + m b gZ b (T] a + rj b ) = Q v . b , 

J| 3 (/3 + I$ 3 ip b + k b ip b = Q Vib , 

m b S a + m b 5 b + dm b 9 a + S b K e ^ b = Qs, b . 

Generalized forces in Lagrangian coordinates are 



(Al) 



(A2) 
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Qe a = a (-yaF a . y + (Z + Z b ) F b . z + (Z + Z a ) F a . z + dF biy ) + 

+ VaFa,Z + S a {F a ,y + F b .y) + (Z + Z b ) F^y + (Z + Z a ) F a .y + 

+M atX - dF btZ + S b F b . y + 9 b Z b F btZ + M btX , 
Qn a = — b — d(p a F b z — 8 a (F ax + F bx ) — (Z + Z b ) F bx + 

777, a 

+ Va ((Z + Z b ) F b . z + (Z + Z a ) F a>z ) + 

+O a (M a , z + M M ) - (Z + + M a + r] b Z b F b , z + M b<y , 

dm b {F atX + <p a F a ,y - r) a F aiZ ) 

Q Va = 1~ dip a F b . y - drj a F b . z + r\ a M b%x - 

m a 

-9 a {M a .y + M b>y ) + r) a M a . x + M a . z + F b . x (d - 9 b Z b ) - r, b Z b F^ y + M M , 
Qs a = —ila {F a , x + F b ^ x ) + 9 a (F as + F b ^ y ) — F a . z — F bjZ , 

Q 8b = Z b (-<p a F b , x + {9 a + 6 b ) F M + F btV ) + tp a M b . y - ri a M b . z + S b F btV + M M , 

Q nb = Z b (-tp a F b .y + (770 + i] b ) F M - F b;X ) - ip a M btX + (6 a + Ob) M biZ - 5 b F b , x + M b . y , 

Q Vb = (r? a + Vb) Mm - + 6b) M b . y + M M , 

Qs b = - (r? a + Vb) F btX + (6 a + 6 b ) F biy - F biI . 

We actually used a further simplified version of these expressions, where we only retain terms of order zero in the 
configuration variables: 

Qe,a - F btV (Z + Z b ) + F a . y (Z + Z a ) + M a>x + M M + y a F UtZ - dF b , z , 
Qr,,a - -F b , x (Z + Z b )-F a , x (Z + Z a ) + M a , y + M b ,y, 

+ dF biX + M ayZ + M M , 
Qs.a — F az F bz: 
Qe,b — Fb,yZb + M byX , 

Qr),b — —F byX Z b + Mb.y, 
Q<p,b — M bjZ , 

Qs,b — —Fbz. 
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